Infrared behavior of the running coupling in scalar field theory 
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We compute the Green function of the massless scalar field theory in the infrared. Applying 
Callan-Symanzyk equation we obtain the exact running coupling for this case by computing the 
beta function. This result is applied using a recently proved mapping theorem between a massless 
scalar field theory and Yang-Mills theory. This beta function gives a running coupling going to zero 
as p 4 in agreement with lattice results presented in Boucaud et al. [JHEP 0304 (2003) 005] and 
showing that the right definition of the running coupling for a Yang-Mills theory in the infrared is 
given in the MOM scheme. 
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I. INTRODUCTION 



' The study of a quantum field theory in a regime where the coupling becomes increasingly large appears to cope 
with a lot of difficulties. The problem relies on the absence of useful methods in this regime. This means that the 
only reliable approach used so far has been lattice computation. 

The model we consider in this paper is a massless scalar field theory. Its relevance relies on the fact that Yang- 
O ■ Mills equations of motion can be reduced to the equation of motion of this theory giving solution for them in the 
infrared. This permits to obtain the Green function for the gluon [l|, Q and the behavior of the beta function using 
Callan-Symanzyk equations. 

For the gluon propagator there is a lot of research activity about, mostly trying to obtain its behavior by solving 
£f~) , quantum Yang- Mills theory on the lattice 0-0] • The scenario that is going to emerge shows that the gluon propagator 
■ goes to a finite value when the momentum goes to zero. On the contrary, the running coupling goes to zero and 
C*") , against the common belief has no fixed point in the infrared. 

The definition used in lattice computation has been given by Alkofer and von Smekal in [3] . This definition is quite 
satisfactory for a qualitative behavior but could be somewhat in disagreement with the real behavior. Indeed, this 
has been seen in a paper by Boucaud et al. 0] that show a different definition for the running coupling using the so 
called MOM scheme and compute its behavior on the lattice obtaining that this goes to zero as p 4 . Alkofer and von 
Smekal expects for the running coupling a fixed point with their definition but this can also goes to zero with the 
expected exponent quite different from the one by Boucaud et al. This opens a question about what should be the 
proper definition for this physical quantity. 

An important answer has been given recently through the analysis of experimental data by a series of papers by 
Prosperi's group These authors take the meson spectrum and get the running coupling by a Bethe-Salpeter 

like formalism. The comparison with analytical methods shows that the experimental data tend to decrease toward 
zero lying well below an expected fixed point in the infrared. 

The current definition of the running coupling, as given in Q, uses the dressing functions of the gluon and ghost 
propagator. One has for the gluon propagator 



for the ghost propagator 



Dip) = ^ (1) 



D G (P) = (2) 



and the running coupling is 

o(p) = J 2 (p)Z(p). (3) 
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We will show consistently that the running coupling to be considered is not this one but rather the one discussed in 
7] with a MOM scheme and going to zero in the infrared as p 4 . 

The paper is structured as follows. In sec. HI] we derive the propagator of the massless scalar field theory and show 
when it solves Yang-Mills equations. In sec. IHII we compare this propagator with lattice results and numerical solution 
of the Dyson-Schwinger equations. In sec. IIVI we see show that our propagator solves Callan-Symanzyk equations 
giving the corresponding beta function. In this way we prove that the running coupling goes to zero in the infrared 
as p A . Finally, in sec. |V] conclusions are given. 



II. TWO-POINT FUNCTIONS FOR MASSLESS SCALAR FIELD AND YANG-MILLS THEORY 



Quantum field theory of a massless scalar field theory is formulated with the partition function 



Z\j} = Jwiexp i J tfxfydtf-ji, 



(4) 



Our approach is to consider the limit A — > oo. This limit has been considered before [TTl - fl9l |. These authors have taken 
^(d<j>) 2 as a perturbation. This choice has the effect to produce a real singular series needing a proper regularization 
with no real chance to get finite results in the given limit. The reason for this relies on the fact that some dynamics 
must be allowed in this limit. This can be obtained by a gradient expansion [l|, 0, [20} . A gradient expansion is 
obtained by rewriting the above functional as 



A 



exp 



i / drx(Vct>y 



(5) 



and considering as a perturbation the gradient part. Now, we note that, in the limit A — > oo, if a dynamics is allowed, 
the term (f> 2 must be of the same order of the term A0 4 . So, we can safely rescale time as t — > VXt and rewrite the 
above functional as 



Z\j}= A#]exp 



exp 



(6) 



and we see that is just a proof of the fact that the limit A — > oo recover a gradient expansion. We need to properly 
rescale the current j to maintain the right ordering but this is just an arbitrary function. The consequence of this 
rescaling of time gives us immediately the conclusion that the strong coupling limit is a semiclassical limit (21.]. So, 
we take <j) = (j) + 8<f> being 



(P + <P =3- 

In the infrared, with the energy going to zero, the following causal approximation does hold [22l . [25 

ft 



being 



#t)« / dt'G(t-t')j(t') 
Jo 



G(t) + G(t) 3 = S(t). 



(7) 



(8) 



(9) 



This is a small time approximation and has been recently recovered in a study of nonlinear waves [24| . It is somewhat 
surprising the proper working of Green functions in such nonlinear systems but this turns out useful to maintain all 
the machinery of quantum field theory. This means that the leading order approximation to hold for our functional 
in the infrared limit is 



Z\j] = exp 



d x\d X2 



X2, 



Sj(x 2 ) 



exp 



d 4 x 1 d 4 x 2 j(xi)A(xi - x 2 )j{x 2 ) 



(10) 



being 



A(zi - x 2 ) - 6 3 {x! - x 2 )[6(ti - t 2 )G(h - t 2 ) + 9{t 2 - h)G(t 2 ~ h)\ 



(11) 
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and 



G(t) = 0(t)/j, [ — I sn 



/j,t, i 



(12) 



being \x an integration constant having the dimension of an energy. So, we showed that a massless scalar theory in 
the infrared take an integrable form. Using the series [29| 



sn[u, i) 



2tt ^ (_i)" e -("+iH 

Kit) ^ 1 + e -(2n+l)n 
v ; n— 



(2n+ 1) 



2#(i) 



being JiT(i) the constant 



a/1 + sin 2 

and after a Fourier transform, we get the full propagator 



1.3111028777. 



G(p) = 



p 2 
n=0 y 



le 



being 



and 



(2n + l) 



7T 2 (-l)»e-(«+2)T 
#2(i) 1 + e -(2n+l) ff ■ 



K(i) 



(13) 



(14) 



(15) 



(16) 



(17) 



From this analysis we see that the theory, in the infrared limit, is trivial having a generating functional as a free 
theory while the propagator is in agreement with the Kallen-Lehman representation |l|, l30( ■ We also note the scaling 
of the propagator with A* that already says to us that the coupling should go to zero as p 4 in a renormalization group 
analysis. We will discuss this in sec lIVI 

This approach is rather general and gives a theoretical framework to treat a quantum held theory in the infrared 
limit. So, an important step is its application to a pure Yang-Mills theory. As already said, there is a lot of activity 
about the solution of this theory in the infrared mostly because there is a serious interpretation problem for the 
spectrum of the light unflavored mesons. For a Yang-Mills theory we will have to solve the equation [31] 



WdpAl - d v {d»Al) + gf abc A b »(d^ - d u A%) + gf abc d»{A b Al) + g 2 f abc f cde A b ^ A d A e , = ft. 



(18) 



being f abc the structure constants of the Lie group and g the coupling constant. When we try to solve for a gradient 
expansion these equations we meet a severe problem. Most solutions are just chaotic [32T - [34j and are useless to build 
a quantum held theory. In order to get a set of solutions to start building a quantum field theory, we can apply a 
mapping theorem recently proved [25l427l | 



Theorem 1 (Mapping) An extremum of the action 



(19) 



is also an extremum of the SU(N) Yang-Mills Lagrangian when one properly chooses A a ^ with some components being 
zero and all others being equal, and A = Ng 2 , being g the coupling constant of the Yang-Mills field, when only time 
dependence is retained. In the most general case the following mapping holds 

A-(x) = r,«<p(x) + 0(l/VNg) (20) 
being constant, that becomes exact for the Lorenz gauge. 



4 



Recent lattice computation in 2+1 dimensions have given strong support to this theorem [28[. So, with this choice, 
the equation for the Green function of a Yang-Mills theory, in a strong coupling limit, reduces to the one of the case 
of the massiess scalar theory with the substitution A —> Ng 2 for a SU(N) group and this is just the 't Hooft coupling 
as should be expected @, H^]- So, finally, e.g. the gluon propagator in the Landau gauge can be written as 

Dfvip) = S ab ~ ^) G(p) + 0(l/VNg) (21) 

being G(p) given in eq.(fT"5"l) with A substituted by Ng 2 . This means that we are approximating the solution of eq.([T8|) 
with 

A%{x) « J d^x'Df^x - x')j h »(x') + 0(l/VNg). (22) 

This latter equation should be checked on lattice computations. 

It would not be difficult to prove that, with this class of solutions of Yang-Mills equations, the ghost propagator is 
the one of a free particle psl]: 

D G (p) = ± + 0(l/VNg). (23) 

giving in our case J(p) = 1 for its dressing function. We can similarly get for the dressing function of the gluon 
propagator, Z(p) = G(p)p 2 . We see that when p goes to zero then we have 

Z(p) a G(0)p 2 (24) 

with G(0) + 1 as can be immediately realized. Then we get the infrared indeces fci = for the ghost and ^2 = 1 for 
the gluon clearly different from those given in 0, [35[ . So one gets for the running coupling as defined by Alkofer and 
von Smekal a(p) — G(0)p 2 that goes to zero in the infrared limit p — > 0. We then get a third index k% = 2 for the 
running coupling. The value of this third index will be the main point of our discussion in sec. IIVI 



III. NUMERICAL RESULTS 



Having obtained an explicit expression for the gluon propagator we have to see how this expression compares with 
respect to numerical results. We will consider two kind of comparison for our aims: firstly we compare it with the 
most advanced lattice results and secondly we will consider a numerical solution of the Dyson- Schwinger equation. 

The problem with Dyson-Schwinger equation is that a scenario has been proposed Q enforcing the view that the 
gluon propagator should go to zero at small momenta, the ghost propagator should go to infinity faster than the 
free particle case and that the running coupling reaches a fixed point in the same limit. Recent data on lattice are 
showing that this view is not correct Besides, data extracted from meson spectrum also showed that the running 

coupling is not going to reach a fixed point bending clearly toward zero 

About numerical solutions of Dyson-Schwinger equation there exist two kind of solutions: on a compact manifold 
[3(1 [33] and D=3+l [38|. For the solution on a compact manifold there are contradictory results but it is clearly 
seen the running coupling bending toward zero. We cannot use these data for comparison. But we can compare with 
the results due to Aguilar and Natale (38|. Solving Dyson-Schwinger equations in D=3+l these authors were able 
to recover the scenario that is presently emerging from lattice computation. They get a gluon propagator going to a 
finite value at lower momenta, the ghost propagator converging to the free one and the running coupling as defined 
by Alkofer and von Smekal going to zero. Then, we must consider these as the proper reference data to compare. 

An important point to be noticed is that there is only one parameter to be computed to compare our propagator 
with numerical data and this is the gluon mass that we write explicitily as m g — (n/2K (i)) fi(N g 2 / 2) i . From this 
value one gets the integration constant /i that is a number to be decided experimentally. 

Firstly, let us see as our propagator compare to the most up to date lattice data We see the results in figHJ At 
lower momenta the agreement is perfect with a gluon mass of 545 MeV. There are volume effects but this are seen 
only in the intermediate energy region. This was also seen in It should be expected that increasing the lattice 
volume should make both curves coincide. 

Finally, we compare our data with the numerical results of Aguilar and Natale [Hj]. In this case we must have 
complete coincidence and this is exactly what happens. This can be seen in figlU The fit is obtained with a gluon 
mass of 738 MeV. From this figure is clearly seen that Aguilar and Natale hit numerically our propagator. With a 
proper fit they could have obtained all the full glueball spectrum. 

The agreements with lattice data and numerical solution are astonishingly good. So, the next step is to draw one 
more important consequence from our solution. 
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FIG. 1: Gluon propagator compared to lattice data in 



IV. CALLAN-SYMANZYK EQUATION AND RUNNING COUPLING 

A proper definition of the running coupling in the infrared is not a trivial matter. An interesting analysis has 
been presented in Q. Lattice computations use the generally accepted definition due to Alkofer and von Smekal [1]. 
It is a obvious matter that a proper definiton can be only obtained having the gluon propagator and solving the 
Callan-Symanzyk equation. This is exactly our aim. Let us point out that in 0] an analysis was carried out through 
lattice computation, assuming a MOM scheme to hold for the running coupling in the form 



c(mom(p) 



4tt 



1 

4tt 



GW(p 2 ,p 2 ,p 2 
(G( 2 V)) 3 



Vg (2 V)) 3/2 



(25) 



being G^ 3 \p 2 ,p 2 ,p 2 ) the three-gluon Green function and G^(p 2 ) the two-point function corresponding to our G(p). 
These authors prove that, on the lattice, cimom(p) ck p 4 . We see that, if G^ 2 \0) — constant, to have the running 
coupling going as p 4 implies G^(p 2 ,p 2 ,p 2 ) oc 1/p and also this is in agreement with Boucaud et al. conclusions. 
Indeed, these authors have given a description of the infrared behavior of a pure Yang-Mills theory in close agreement 
with the one that is presently emerging from lattice computations (see e.g. [42| and refs. therein). 

The propagator of the massless scalar theory we have derived above must satisfy the Callan-Symanzik equation 



dG(x,t) n „ s dG(x,t) s~i/ \ n 
V — tt^ 2 - + — + 2 7G(x, t) = 



and setting x = we have immediately 



9G(t) n ,^dG(t) n „, . n 



or, working with momentum, one has 



We see immediately that this is true if 



dG(p) 
<9/i 



/3(A) 



dG{p) 



2 7 G(p) = 0. 



(26) 



(27) 



(28) 



/3(A) = 4A 



(29) 
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FIG. 2: Gluon propagator compared to numerical data in [38 



and 

7 = -1 (30) 

givi ng us the behavior of the running coupling we looked for. This result is in agreement with recent analysis 
(39L |40j where such a form of beta function in the strong coupling limit was postulated. Similarly, from AdS/CFT 
correspondence a similar conclusion was drawn (4lj . 
We know that 

0(A) = M £ (3D 



giving immediately 



A(p) = A 4- (32) 



Now, with the substitution A — > Ng 2 that maps the scalar theory on the Yang-Mills theory we get immediately the 
result given in [7|. This permits us to draw immediately the conclusion that the definition cxmom{p 2 ) is the proper 
one for the running coupling in the infrared. This is expected to go to zero with a fourth power law reaching no 
fixed point. This gives index k% being 4. Such a result can be understood if we recall what we have said above about 
the ghost propagator. The ghost in the infrared limit decouples from the gluon and then, a definition of the running 
coupling implying the ghost propagator has no physical meaning. 
Similarly, we have 

1 M dZ 

^2ZdJl (33) 

being Z the renormalization constant of the field. We see immediately that Z = p 2 //j? being this the expected scaling 
for the field. We just note that in the infrared limit particles are not gluons but rather the particles in the mass 
spectrum of the theory that should be properly called glueballs. 
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V. CONCLUSIONS 



The conclusions to be drawn describe a scenario completely different from the one discussed in literature in these 
latter years. This is due mostly to lattice data that have shacked during these years the common beliefs that were 
going to form on known theoretical methods. 

Presently, evidence is mounting that the gluon propagator reaches a finite value for momenta going to zero. In this 
paper we have shown that the proper definition of the running coupling in the infrared should be taken in a MOM 
scheme. This is shown to go to zero, not to a fixed point, with a fourth power of momentum. Finally, this result 
agrees well with the view that the ghost in the infrared decouples from the gluon and behaves as a free particle. 

As this scenario is emerging from numerical solutions of the Yang-Mills quantum field theory and from phcnomcno- 
logical analysis, increasing confirmations in the future years have to be expected. 
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